Medium-induced parton splitting kernels from Soft Collinear Effective Theory 

with Glauber gluons 
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We derive the splitting kernels for partons produced in large Q 2 scattering processes that subse- 
quently traverse a region of strongly-interacting matter using a recently-developed effective theory 
SCETg. We include all corrections beyond the small- a; approximation, consistent with the power 
counting of SCETg- We demonstrate how medium recoil, geometry and expansion scenarios, and 
phase space cuts can be implemented numerically for phenomenological applications. For the simpli- 
fied case of infinite transverse momentum kinematics and a uniform medium, we provide closed-form 
analytic results that can be used to validate the numerical simulations. 



I. INTRODUCTION 

The suppression in the production rate of energetic 
leading particles and particle correlations due to final- 
state interactions in reactions with ultra-relativistic nu- 
clei is among the best-known experimental discoveries at 
the Relativistic Heavy Ion Collider (RHIC) 0-|§, and 
now at the Large Hadron Collider (LHC) §-0! This jet 
quenching phenomenon also provides one of the strongest 
pieces of evidence for the creation of dense strongly- 
interacting matter in such collisions . Recent advances 
in understanding the modification of partons and par- 
ton showers in Quantum Chromodynamics (QCD) media 
have come from the inclusion of jets in the theoretical Q- 
ijlfjf and experimental []l7j-[E3) analyses. Jet observables 
are more sensitive to the underlying theoretical model as- 
sumptions and the properties of the QCD medium when 
compared to leading particle measurements |24| . 

In recent years effective field theories (EFT) have be- 
come a powerful modern tool for jet physics. In particu- 
lar, Soft Collinear Effective Theory f|j-|| (SCET) is an 
effective theory for QCD that describes the dynamics of 
highly energetic partons. It has been successfully applied 
to improve the theoretical accuracy in the evaluation of 
high energy cross sections at lepton [p9|-|3~il and hadron 
colliders [p3|- 35 



The first step in constructing an effective theory for 
jets propagating in a QCD medium was done in Ref. [ [36| , 
where the SCET Lagrangian was extended by adding a 
term that describes the interaction of a quark jet with 
gluons that have momentum purely transverse to it, tra- 
ditionally referred to in the literature as Glauber gluons. 
As an application of effective Lagrangian derived in |36| , 
the probability density of quark jet broad- 
39(| was re-derived as an expectation value 
of Wilson lines, which later the authors evaluate using 
AdS/CFT correspondence. In Ref. J5} the Yang-Mills 
part of the collinear SCET Lagrangian was coupled to 
Glauber gluons, which allowed to perform calculations 
for parton splitting processes in the medium. The re- 
sulting effective theory was called SCETg, where "G" 
stands for Glauber gluons. In that paper a detailed con- 
nection was made between calculations in SCETg and 
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the evaluation of jet broadening j38[ p9f and medium- 
induced quark energy loss in the Gyulassy-Levai- Vitev 
approach |38[ III]]. The gauge invariance of the physics 
results was explicitly demonstrated for three different 
gauge choices. One of the medium-induce d sp littings, 
namely q — >• qg has been calculated in Ref. |4(| beyond 
the soft emitted gluon approximation. There are three 
additional splittings: g — > gg, g qq and q — > gq. To 
complete the derivation of all medium-induced branch- 
ing processes without the assumption of a soft final-state 
parton and to understand the correction that arise from 
the finite parton scattering kinematics, branching kine- 
matics, and recoil of the constituents of the QCD medium 
is the main goal of this Letter. 

The rest of this Letter is organized as follows: in 
section || we discuss the theoretical framework SCETg 
for our calculation and demonstrate how the vacuum 
Altarelli-Parisi splitting kernels can be derived in Soft 
Collinear Effective Theory. The derivation of the full 
splitting kernels for in-medium jet production with final- 
state interactions is discussed in section III . We elucidate 



the relation to early soft gluon approximation results 
and provide analytic formulas for simplified kinematics 
and medium geometry scenarios. Numerical control over 
the newly-derived medium-induced splitting intensities is 
demonstrated in section IV . In this section we also quan- 



tify the effects of large- a; corrections, finite kinematics, 
and medium recoil. A brief summary and outlook is pre- 
sented in section |v|. 



II. THEORETICAL FRAMEWORK 



An effective theory, well-suited to describing the prop- 
agation of jets in the medium, has been recently devel- 
oped in Ref. j!o| . The Lagrangian of this EFT is given 
by the sum of the SCET Lagrangian p5|-p8| and a term 
that specifies the interactions of collinear partons in QCD 



2 



matter: 



-£SCET G (£n , A n , A G ) — ^SCET (£nj A n ) + Cq (£n, A ni A G ) 



= fl, 0, 0, 1) , = (1, 0, 0, -1). The action in Eq. 
Eq. <H) is given by Lagrangian of SCETg: 



S = 



d x£sCET G 



(9) 



(1) 



In Ref. ||FJ the vertexes r^^r^^ have been de- 
rived for three types of gauge-fixing conditions: covari- 
ant, light-cone and hybrid gauges. In the first case we 
gauge-fix both the physical collinear gluons as well as 
the Glauber gluons in the covariant gauge. The second 
choice corresponds to gauge-fixing both fields using the 
light-cone gauge. The third choice, which appears to be 
the most convenient from the practical point of view, cor- 
responds to a light-cone gauge for collinear gluons and a 
covariant gauge for the Glauber gluons. This is a legit- 
imate choice from effective theory point of view, since 
we are allowed to gauge-fix separate gauge sectors inde- 
pendently. Another way of justifying this gauge choice is 
factorization between the splitting and the elastic scat- 
tering. In this hybrid case both the collinear Wilson line 
W and the transverse gauge link T [ff2|-[44| vanish. Gauge 
invariance of the physics results for the in-medium elas- 
tic scattering and radiative energy loss was demonstrated 
in jio) , providing a cross-check on the approach and the 
newly-derived Feynman rules. It is interesting to note 
that the same effective theory SCETg is relevant for de- 
scribing the Drell-Yan process, as shown in Ref. p5[ . 

We start from amplitudes for the parton splitting pro- 
cesses: 



The Lagrangian of SCETg §| @ is S iven in Ec l- © and 
it evolves the created jet and describes the parton split- 
ting processes and the interaction of the parton shower in 
the medium. The amplitude with qq O qq is not shown 
explicitly. 

Restricting ourselves to the SCET Lagrangian with- 
out Glauber gluons, we first verify that at tree level we 
recover the Altarelli-Parisi splitting kernels [^8| , which 
have been originally calculated in full QCD: 



dN 
dx G? 2 kj_ 

dN 
dxd 2 h± 
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dx d 2 k± 
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dx d 2 k± 



a s I + (I — a;) 2 1 
Of 
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9^99 



2tt 



(10) 
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2tt 2 \ x 1 — x 



1 



(ii) 



^T R (^ + (1-^)^,(12) 



2tT 2 " 



dN 

dx G? 2 k^ 



(x 1 — x). 



g^gg 



(13) 



A q ^ qg = (q(p)g(k)\Te tb Xn(xo)\q(Po)), (2) 

A g ^ gg = (g(p)g(k)\Te* s B Xc (x )\g(p )), (3) 

V>, g = (q(p)q(k)\Te iS B Xc (x )\g(p )) , (4) 

A q ^g q = (g(p)q(k)\ Te lS Xn(x ) \q(Po)) , (5) that we derive in section fn 



We note that we are interested in real splitting pro- 
cesses away from the singular end points x — and 
x = 1. In all expressions above the transverse mo- 
mentum k± and the lightcone momentum fraction x = 
k + /Po = k + /(p + + k + ) are for the second final-state 
parton. The parent parton has no net transverse momen- 
tum and = —p±- Note that Eq. (10) and Eq. (13) 
are interchangeable under x — > 1 — x, whereas Eq. (11) 
and Eq. (^|) are symmetric under this substitution. The 
same symmetries hold for the medium-induced splittings 



where B are collinear gauge invariant SCET fields |4q , 
f47{ and the momentum four-vectors, such as po = p + fc, 
are parametrized in the standard way, consistent with en- 
ergy momentum conservation and the on-shell condition 
p 2 = k 2 = 0: 



Po 
p 
k 



Po 



x(l- x)p£ 



(1 - x)p+ 



(1 - x)p^ 



k 2 

x Pq J 4- i K i 



(6) 
(7) 
(8) 



We use square brackets to indicate the light-cone nota- 
tion, which we define for arbitrary four-vector q in the 
following way: q = [q + ,q~,q±] = [n-q,n-q, q ± ] and 



In this paper we use the following terminology: the 
double differential distribution dN/dxd 2 ]i± we call a 
splitting kernel, xdN/dx we call a splitting intensity and 
dN I dx we call a differential emitted parton number dis- 
tribution. This terminology applies to both vacuum and 
medium-induced splittings. The x— dependent part of 
the vacuum splitting kernel we call a splitting function. 
Since the medium-induced kernel has a more complicated 
k± , x correlation structure compared to the simple fac- 
torized form in Eq. (10) - Eq. (13) we avoid definition 
of a similar term in the medium. 



III. MEDIUM-INDUCED PARTON 
SPLITTINGS 

To describe the collisional and radiative processes 
for partons propagating in QCD matter, both single 
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dN 
dx 



+ 2Re 




+ /- 



x i- 



FIG. 1: Feynman diagrams contributing to medium-induced 
splittings at first order in opacity. Red lines corresponds to 
Glauber gluons. The kinematics and topology are common to 
all splitting processes: q — » qg, g — > gg, g — > qq, q — > gq. 



and double Glauber gluon exchanges between the jets 
and the constituents of the medium must be consid- 
ered [j|o| [II . 50 . The calculation to first order in opacity, 
which takes into account the contribution from the split- 
ting induced by the interactions along the trajectory of 
the parent parton and the dominant interference with 



trated in figure |[ We do not specify the parent and 
daughter parton flavors since the topology and kinemat- 
ics are the same for the splitting processes enumerated 
in Eqs. (|l(]) - ([l3]). Consequently, all results can be ex- 
pressed in terms of universal transverse momentum vec- 
tors A±, B±, Cj_, D± and interference phases Hi, Og, 
defined in Gol: 



A± 


= fcx, 


B± = 


D± 


= k ± 








B 


fix 


-o 2 = 


p+x(l 







(14) 



n 2 - o, = 



C 2 — B 2 , 



Pqx(1 - x) 



n 4 = 



Pqx(1 - x)' 



Pqx(1 - x) ' 



n« = 



P$x(l - x) 



(15) 



the splitting induced by the large Q process, is illus- 



where p^ = p + + k + and the parent parton has no net 
transverse momentum. 

For completeness, we first present below the result for 
the q — > qg splitting, calculated in Q and shown to be 
gauge invariant: 



dN 



dxd 2 k± 



a s l+(l-x) 2 
Of 



2tt 



dAz 



1 dcr " lcdlum 

d V 4 

(J ei a qx 



B 2 



B_ 

B 2 



g^gg 

: (1 - cosKfir - n 2 )A*]) + (2^ - ^ - fr) ( X - cos ^i - ^)Az]) 
-|r-§r (! - «»[(fi a - fi 3 )A*]) + ^ - ^) (1 - cos[fi 4 Az]) 

- cos[0 5 Az]) + - |r) (1 - oobKOi - naJA*]) 



(16) 



where A g (z) is the scattering length of a gluon in the 
medium and (l/cr e ;) d<7 e ™ edium /d 2 q_|_ stands for normal- 
ized elastic scattering cross section of a parton in the 
medium. Even though this quantity varies when parton 
is a quark or a gluon, in the high energy limit, when the 
t— channel dominates the elastic scattering, this normal- 
ized cross section does not change significantly. 

Using the Feynman rules of SCETq in the hybrid 
gauge and the Feynman diagrams exactly analogous to 



the case of q — > qg splitting considered in pQ ] and shown 
in figure we derive the remaining parton splittings in 
the medium. The calculations are non-trivial and facili- 
tated by intermediate results in . As discussed in sec- 
tion |l|, the medium-induced splitting for q — > gq can be 
obtained from Eq. ( |l6| ) with the substitution x — > 1 — x. 
Here, we skip the explicit expression for brevity. The re- 
maining two splittings from a parent gluon are as follows: 
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where X q (z) is the scattering length of a quark in the 
medium and same comment applies to the quantity 
(1/o-ej) da£ edium /d 2 <lx as after Eq. ©. Note that up 
to the overall vacuum-like splitting functions and color 
factors reflected both in the mean free paths (quark ver- 
sus gluon) and the corrections relevant beyond the small- 
x approximation, the structure of the answers is very 
similar. The symmetry of g — ¥ gg,g — > qq splitting ker- 
nels under x — > 1 — x is most easily verified explicitly by 
realizing that the parton scattering cross section in the 
medium is invariant under — > —q_i_. 

The basic features of the medium- induced parton split- 
ting kernels are: 

• In QCD, for parent quark they factorize from the 
hard scattering cross section up to a standard inte- 
gral convolution pOfl . For parent gluons non-trivial 
spin correlation are present analogous to the vac- 
uum case 1491. 



• They are proportional to their vacuum Altarelli- 
Parisi splitting functions Q . 



The in-medium splittings are gauge-invariant, as 
they should be, since the underlying jet production 
process itself is gauge-invariant 



• The splitting kernels depend on the properties of 
the QCD matter and vanish when the size or den- 
sity of the medium vanish. The functions derived 
here are only valid for final-state interactions 

It is instructive to verify that in the small- a; limit only 
two of the four medium-induced splitting intensities sur- 
vive and this allows for the standard energy loss inter- 



pretation of jet quenching: 




0{x)] 
0(x)} 
0(x 2 )} 
0(x 2 )] 



1 dcr ™ edium 



2k± ■ q-j 



(fcx - Qx) 2 A ' 
1 — cos r Az 

x Po 



(18) 



In this limit the interference structure for all medium- 
induced splitting intensities is the same. Furthermore, in 
the small- a; limit the last two splittings are suppressed 
(0(x)) relative to the first two. We keep the first correc- 
tion for numerical comparison only. The color structure 
for the in-medium interactions also simplifies in this limit 
and is determined by the flavor of the small- a; parton in 
the final state. Specifically, the first two in-medium split- 
tings are proportional to 1/A g and the second two are 
proportional to 1/A ? . In deriving these results, we have 
used relation: X q /X g — Ca/Cf, which follows from the 
leading order perturbation theory approximation. As ex- 
pected, in the small- a; emission limit our results coincide 
exactly with the intensity derived (or neglected when the 
leading term is 0) in [E{J . 

In section IV we will study numerically the in-medium 
splittings derived here with an emphasis on going beyond 
the traditional small- a; approximation and on including 
medium recoil. The remaining part of the current section 
is devoted to deriving analytic formulas for the inclusive 
splitting intensity x(dN/dx) under certain idealized as- 
sumptions. This will, in turn, allow us to obtain fully 
analytic formulas that can be used to benchmark the re- 
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alistic numerical calculation. 

A useful starting point for integrating the splitting ker- 
nels analytically over the transverse momenta is the fol- 
lowing master formula: 



d 2 k ± d 2 q 



1 da. 



2k, 



°ei cPq_L k 2 ± (k± - q ± ) 2 
[l-cos(a(k ± -q ± ) 2 )] =f[a»% (19) 



where a specific form for elastic cross section was assumed 
as explained below and shown in Eq. (^l]), and function 
f[x] equals: 



2tt 



7e + In (a; 
sin(a;)Si(x) 



■ sin(a;) — cos(x)Ci(ir) 



(20) 



Two assumptions are already made at this level. First, 
we took the limits of integration on fej_, to infinity. In 
reality, phase space cuts affect the cross section and we 
study this effect numerically in the next section. Second, 
we neglected the recoil effect in the medium. In that 
approximation the normalized cross section equals: 



1 da P 



a e i d 2 q±_ n{q\ + /i 2 ) 



(21) 



The effects of finite medium recoil are also studied in 
the next section numerically. It turns out that using 
Eq. (19) it is possible to calculate k±,q ± integrals in 
all in- medium splittings Eq. (|l6|)-Eq. (17|). The result 
is rather compact and can be expressed in terms of the 
function f[x] defined in Eq. (20): 



dN 
dx 

dN 
dx 

dN 
dx 

dN 

d,x 



q^qq 



q^qq 



a s „ l + (l-x) 2 f dAz M + ~ xf] ~ 



2tt 2 



X 

x 



x 



A fl (z) 
1-x 



+ x(l - x) 



2ir 2 

dN\ 
dx J 



T R (x 2 + (1 - x) 2 ) 



q-^qg 



dAz 



(x — > 1 — x), 



N? 



dAz f[px 2 } +fW -x) 2 ]+f[P] 
\ g (z) 2 



N 2 1 
" - [f[px 2 ] + f[p(l - x) 2 ]) - j^—f\p\ 



(22) 
(23) 
(24) 
(25) 



r 



where the superscript oo stands for infinite limits of in- 
tegrations for k±,q ± , and: 



\i 2 Az 



Pox(l - x) 



(26) 



In order to perform the remaining Az integral, one 
has to specify the geometry and an expansion scenario 
for the QCD medium that the partons traverse. Even 
the simplest realistic model of the medium, which in- 
cludes the Glauber nuclear geometry and Bjorken ex- 



pansion, requires numerical evaluation of the in-medium 
splitting kernels. One can validate the numerical simula- 
tion techniques by comparing the results to closed-form 
analytic formulas for uniform QCD matter, where the 
parton mean free paths X q , X g and trajectory length L 
are fixed. Thus, it is instructive to have analytical for- 
mulas for splitting intensity for uniform QCD matter. By 
assuming constant \ qtg as a function of Az, and integrat- 
ing expressions in Eq. (p2|)-Eq. (p5|) over < Az < L, we 
obtain the following differential splitting intensities: 



dN_ 

dx 

dN 
dx 

dN 



q->-qg 

oo, static 



g^-gg 

oo, static 



a s 1 + (1 - x) 2 L , , 
x—C F — ( g[j] + g[-y{l - x) 1 

Z7T X An 



1 

N 2 



a, „ / x 



+ + x(l - x)) — (g[ 1X 2 } + <?[ 7 (1 - x) 2 ] + g[i\) , 

X I A n 



2N 2 2 

(ahx 2 ] + gh(i - x) 2 ]) - jfT^gb] 
i 



(27) 
(28) 
(29) 



where the function g is given by: 

7T 

w g[w] = — (1 — cos(u>)) + (te — l)w + w ln(w) 



and 7 is defined as: 



fi 2 L 



+Si(w) cos(w) — Ci(w) sin(iu), 



(30) 



Pqx(1 - x) 



(31) 
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FIG. 2: The intensity spectrum x{dN/dx) for infinite phase space cuts and neglec ting nuclear recoil is shown as a function of 
the splitting parameter x. Comparison of the analytic formulas in Eq. @-Eq. (||) (solid lines) to a numerical integration 
method (dashed lines) are presented. We also illustrate the difference between the full in-medium splitting results and the 
small- a; approximation on the example of a parton of initial energy Eq = lOOGeV. The medium parameters are set to: 
fi = 0.75 GeV, A g = 1 fm, L = 5 fm for definiteness and the scattering length is independent of Az. 



The intensity spectrum for the last splitting q — > gq can 
be obtained from substitution x — > 1 — x in the q — > qg 
splitting and is given in Eq. (|25|). 



IV. NUMERICAL RESULTS 

In this section we study the effects of kinematic cuts 
and recoil of the medium by evaluating dN/dx numer- 
ically. In so doing, we demonstrate control over the 
numerical evaluation, keeping in mind that future ap- 
plications will require such approach to incorporate the 
finite kinematics, the spatially non- uniform and time- 
dependent density of the QCD matter, and the recoil 
of the in-medium partons. For each splitting we con- 
sider the full result given by Eq. (16) - Eq. ( |i~7| ) and 
compare it to the small- a; limit presented in Eq. ( |l8| ). 
In this paper we consider a medium of uniform density 
for simplicity and set the parameters of the simulation 
as follows: the typical inverse range of the parton scat- 
tering in the medium is fj, = 0.75 GeV, the size of the 
QCD medium is L = 5 fm, the gluon mean free path in 



matter is X g = 1 fm, and the parent parton energy is 
E Q =p+/2= 100 GeV. 



For infinite limits of the k± 7 q ± integrations, ignor- 
ing the medium recoil effects and assuming static QCD 
matter, we checked numerically our analytic formulas in 
Eq. (^) - Eq. (p9|). We found perfect agreement that 
validates the numerical integration methods and the an- 
alytic results. This can be seen from figure 0. Solid lines 
represents the analytic results of Eq. (pTj) - Eq. (|2S|) . 
Dashed lines represent numerical results. Our conclu- 
sions are valid for both the full in-medium splitting in- 
tensity x(dN/dx) and its small- x limit. Note that for 
such comparison to be possible we have retained the sub- 
leading 0(x) term for the g — > qq and q — > gq processes. 
As expected, the deviation between the full in-medium 
splittings (red and blue lines) and their small-x approx- 
imation (green and black lines) is the largest as x — > 1. 
For intermediate x ~ 0.5 the deviation is on the order of 
a factor of 2 and changes sign. 

In figure [| we present the comparison of the splitting 
intensities without transverse momentum cuts and with- 
out parton recoil in the medium (solid red curve) to three 
different cut scenarios. In all thr ee cases we use the 
same cut on fc max = y/ Q 2 x(l — x), which is unambigu- 
ous, and we choose Q = Eq. The three scenarios for the 
cut are: o) the dashed green curve corresponds to 
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No cuts, no recoil 

^ma^^O^ 2 ) 172 - k max =E oW 1 - X )) 1/2 ' no recoil 

q rTiax =(2E Q u.) 1/: , k ma< =E (x(1-X)) 1/2 , no recoil 
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Full kinematics cut, recoil m =1 GeV 
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Full kinematics cut, recoil m =1 GeV 
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FIG. 3: Illustration of the effect of phase space cuts and medium recoil on the medium-induced parton splitting. The same 
QCD medium parameters and initial jet energy as in figure H are used. 



<7max = b) the dashed blue curve corresponds 

to g m ax = \/2^Eq, c) the dot-dashed black curve corre- 
sponds to the exact phase space, given by 2 —> 2 scatter- 
ing available phase space. Finally, the solid black curve 
includes the recoil effect which is calculated by substi- 
tuting the normalized cross section in Eq. (El]) by the 
2^2 i-channel differential cross section, which can be 
found in Eqs. (3.2-3.3) in Ref. @. From the definition 
of cut on fcj, for small x we have fc max ~ Q\fx 0. 
From figure [| one can see that for small x the cut on k± 
is the only one that affects the splitting intensity, since 
all three types of cuts on q ± give practically identical 
results. As far as the intermediate x region is concerned, 
the cut on k± docs not play a significant role since from 
the definition fc max (a; ~ 1/2) ~ Q/2, thus the observ- 
able difference must be attributed to the cut on q ± for 
this region. The kinematic cut on however, can lead 
to a factor of 2 variation of the in-medium parton split- 
ting intensities at intermediate x. Note that for interme- 
diate x the third cut on q±, which corresponds to full 
kinematics but retains the l/tfj dependence of the scat- 
tering cross section Eq. (plj), agrees perfectly with the 
uncut solid red curve. In this case, cuts alone (in the 
sense of full kinematics) do not affect significantly the 
in-medium branching processes. We find that what af- 
fects the splitting intensity is the deviation between the 
exact scattering cross section from Q and the approx- 



imate power-law form in Eq. (|2l|). This is illustrated 
in figure || by the solid black curve that pushes the in- 
tensity of the medium-induced branching processes down 
when compared to the dot-dashed black curve. We finally 
note that if one wishes to simplify the calculation and 
use the approximate form Eq. IX2J\) for in-medium parton 
scattering the most adequate transverse momentum cut 
would be Qmax = vJiEq. Finally, for x — > 1, we find that 
fcmax ~ Qvl — x 0, and all the splitting intensities 
with phase space cuts turn over at large enough x, which 
is not visible in figure || because we do not plot values of 
x very close to 1. 

Numerically, all effects: finite x, phase space cuts, re- 
coil effect are of the same order at high energies. In ad- 
dition, we observe that parton recoil, similar to finite x 
corrections appears at intermediate x, while phase space 
cuts play role both for small x and intermediate x. 

V. CONCLUSIONS 

In this Letter we derived the medium-induced parton 
splittings for quarks and gluons produced in large Q 2 
scattering processes that subsequently traverse a region 
of dense QCD matter and undergo final-state interactions 
using a recently constructed effective theory SCETq [ fl0| . 
Our results include both the contributions from the in- 
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medium parton scattering and their dominant interfer- 
ence with the vacuum-like branching processes 0, [50| . 
These formulas are valid for an arbitrary splitting pa- 
rameter x and include all leading terms consistent with 
the power counting of SCET G [M. Our results for the 
q — > qg, q — > gq, g — > gg, g — > qq, including Landau- 
Pomeranchuck-Migdal coherence and interference effects, 
are presented in Eq. (|l^)-Eq. ( |l7| ) and are the main new 
findings reported in this Letter. We verified in Eq. ( |l8| ) 
that in the small- a; approximation our formulas simplify 
considerably and reduce to the soft gluon emission results 
of the Gyulassy-Levai-Vitev approach to inelastic parton 
interactions in dense QCD matter |4l], |50| . In this case, 
to leading power in 1/x, only two medium-induced split- 
ting kernels survive that do not change the flavor of the 
leading parton and have a natural interpretation in terms 
of parton energy loss. 

Neglecting the recoil of the partons in the QCD 
medium and ignoring the phase space cuts, we derived 
fully analytic formulas for the in-medium splitting inten- 
sity, given in Eq. (p7|)-Eq. (|29|). We see the main util- 
ity of these formulas as a convenient cross check for the 
numerical simulations, which we have demonstrated in 



this Letter. Our full results provide a basis for further 
improvements in the jet quenching phenomenology |]l0|| - 
p"5| by including the following effects: a) finite x correc- 
tions in the q — > qg and g — > gg splittings, consistent 
with the power counting of SCETq , b) new splittings for 
g — > qq, q — > gq from coherent final-state interactions, c) 
exact parton recoil kinematics and, d) exact phase space 
cuts. All of the above effects can be of the same order, 
as one can see from our numerical simulation results in 



section [V. This program of improving the theoretical 
accuracy of jet quenching simulations is especially inter- 
esting in light of recent RHIC and LHC heavy ion results 

[@"@> anc ^ we P^ an *° snow nrs t phenomenolog- 

ical applications in the near future. 
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